Bright entangled light from two-mode cascade laser 
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We show that a two-mode three-level cascade laser driven by external coherent fields generate 
intense entangled light. It turns out that external driving fields which are at resonance with the 
cavity modes substantially improves the intensity of the two-mode light in the cavity in a region 
where the squeezing and entanglement is significant making the system under consideration a viable 
source of bright squeezed as well as entangled light. 
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I. INTRODUCTION 

Generation of entanglement has recently attracted 
great interest as it plays a key role in quantum infor- 
mation processing |l|-|5(. Particularly, much attention 
has been paid to generation of continuous-variable en- 
tanglement as it might be easier to manipulate than the 
discrete counterparts, quantum bits, in order to perform 
quantum information processing. In general, degree of 
entanglement degrades as it interacts with the environ- 
ment. On the other hand, the efficiency of quantum in- 
formation processing highly depends on the degree of en- 
tanglement. As a consequence, it is desirable to generate 
strongly entangled continuous-variable states which can 
survive from the environmental noise. 

Schemes for generating continuous- variable states have 
been realized in optical parametric oscillators |6j,|7[. More 
recently, two-mode three level cascade lasers have been 
proved to be a source of macroscopic bipartite entangled 
states. For example, Xiong et. aZ.:8j, have demonstrated 
macroscopic entanglement in a driven two-mode three 
level laser when the atoms are injected at the lower level, 
applying the entanglement measure proposed by Duan 
et al.^. Tan et a/.[l(| extended this work and stud- 
ied the generation and evolution of entanglement using 
the Wigner representation. The steady state entangle- 
ment in a two-mode three-level laser, where the atomic 
coherence is induced by initially preparing atoms in co- 
herent superposition of the top and bottom levels has 
also been studied [Tl|. The effects of coupling this sys- 
tem with optical parametric oscillator [l^] and two-mode 
squeezed vacuum reservoir |13[ on the steady state two- 
mode squeezing as well as entanglement has thoroughly 
been studied more recently. It is found that the entangle- 
ment in such systems is directly related to the two-mode 
squeezing. 

In this paper, we propose a scheme to produce bright 
entangled light from a two-mode three-level laser driven 
by two coherent fields. Three-level lasers are long known 
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for a source of squeezed light where the crucial role is 
played by atomic coherence which can be induced by 
initial preparing the atoms in a coherent superposition 
of the top and bottom levels [l4- 19 j or coupling these 
levels by external driving fields On the other 

hand, the effect of driving the cavity mode by external 
coherent light has been studied recently [l9j]. In the ex- 
perimental setting, driving the OPO cavity with seed 
waves in a way to improve the intensity of the signal 
and idler modes which resulted from the down conver- 
sion process has been realized [23| . In this work we show 
that the driving coherent fields considerably enhance the 
mean photon number, without affecting the degree of en- 
tanglement obtained from the laser system, in a region 
where the squeezing and entanglement is significant. This 
makes the system under consideration a source of intense 
squeezed as well as entangled light. 

We derive the pertinent master equation in the lin- 
ear and adiabatic approximation schemes. The resulting 
master equation is used to obtain equation of evolution 
for the first- and second- order moments of the cavity 
mode operators. Using the steady state solutions of these 
equations we study the two-mode steady state squeezing 
in sum and difference fields. Moreover, applying the same 
solutions, we analyze the entanglement properties apply- 
ing the entanglement measure proposed Duan et al. [9J. 
Finally, we also calculate the mean photon number of the 
two-mode light. 



II. HAMILTONIAN AND MASTER EQUATION 

We consider nondegenerate three-level cascade atoms 
injected in to a cavity coupled to a vacuum reservoir. 
The atoms are initially prepared in a coherent superpo- 
sition of the top and bottom levels in order to induce 
atomic coherence to the system. The atoms are injected 
into the cavity at some constant rate r a and removed af- 
ter they decay spontaneously other than the middle and 
intermediate levels. We assume that, the transition from 
the upper energy level \a) to the intermediate level \b) 
and from level \b) to the lower energy level |c) are taken 
to be resonant with the cavity modes, whereas the tran- 
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sition | a) —> \c) is dipole forbidden. In turn the cavity 
modes are driven by two external coherent fields having 
the same frequency as the respective cavity modes. 
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FIG. 1: (a) Scheme of a two-mode three-level cascade laser 
driven by two coherent fields having the same amplitude, e 
but different frequencies, ui a and uib- The gain medium is 
ensembles of three-level cascade atoms, (b) Three-level atom 
in a cascade configuration. 

The interaction of the cavity modes with the external 
driving fields and with a single three-level cascade atom 
is described, in the rotating wave approximation and in 
the interaction picture, by the Hamiltonian 

H = ie{a) - a + S f - b) 

+ ig@\b) (a\ + &t| c ) (b\ - \a) (b\a - \b) (c\b), (1) 

where e is the amplitude of the external driving field as- 
sumed to be the same for both fields, g is the atom-cavity 
mode coupling constant assumed to be the same for both 
transitions, a and b are the annihilation operators for the 
two cavity modes. In this work, we take the initial state 
of a single three-level atom to be 
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(2) 



IM0)> =C a \a 
and the corresponding density operator is 

/U(0) = p^\a)(a\+p^\a)(c\+p^*\c)(a\+p^\c)(c\, 

(3) 

where pi? = |C Q | 2 and pi? = |C C | 2 are respectively the 
probabilities for the atom to be initially in the upper and 
lower levels and pi? = C a C c = pi?* represents the initial 
atomic coherence of the atom. 

We are interested in the evolution of the cavity modes 
only. This can be achieved by tracing the atom-cavity 
modes density operator over atomic variables. The time 
evolution of the density operator for the cavity modes 
pit) can, in general, be written as 



-p(t) = -iTr A [H,p AR (t)}. 



(4) 



Now employing Eqs. ©-(Ill) and taking into the damping 
of the cavity modes by the vacuum reservoir, we obtain 
the master equation for the cavity modes, in the linear 



and adiabatic approximation schemes, to be of the form 
d 



Jt P = e^p 



pa" — ap + pa) 
+ e(& t p - -bp + pi) 



(2a) pa — aa) p — pact') 



~(2apa) 



a'ap — pa) a) 



Ap?J 



{ptfP +a t S t p- 2a t p6 t ) 



(pab + abp — 2bpa), 



where 



A = 2.g 2 r a / 7 2 



(5a) 



(5b) 



is the linear gain coefficient, 7 is the spontaneous decay 
rate assumed to be the same for all the three levels, k is 
the cavity mode damping constant which we assumed it 
to be the same for each cavity mode for convenience. 



III. SQUEEZING IN THE SUM AND 
DIFFERENCE FIELDS 



In this section, the squeezing properties of the two- 
mode light in the cavity produced by the two-mode three- 
level cascade laser is analyzed. The squeezing properties 
a two-mode light can be investigated by introducing sum 
and difference fields of the two cavity modes and analyz- 
ing their respective quadrature variances. 

We define the sum operator of the two modes as 

c=-L(a + S). (6) 
The corresponding quadrature operators are 



c+ = c 1 + c, 



i(c f 



(7) 



Based on these definitions of the quadrature operators 
the two-mode light represented by the sum operator is 
said to be in a two-mode squeezed state provided that 
either of the variances of these operators, Ac 2 ^ or Acl, 
should be less than that of the vacuum level, which is 
unity. The variances of the quadrature operators can be 
expressed as 



A4 



(4) 



{c±) 



(8) 



We wish to calculate the quadrature variances at steady 
state. Thus, using the steady state solutions obtained at 
the Appendix, the steady state quadrature variances are 
found to be 



Ac 2 = 



A 2 (l - if) + (2k + A + At])(2k + At]± Ay/1 - if) 



2(2k + At]){k + Arf) 



(9) 
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FIG. 2: Plot of the quadrature variances Ac?, and Ad+ versus 
k for Tf = 0.1 and A=100. 

Similarly we can define the difference operator of the 
two modes as 

d=-L(&-6) (io) 

with corresponding quadrature operators 

d+ = $ + d, <L = *(# - d). (11) 
The variances of these quadrature operators are: 

a4 = (4) - (d ± f. (12) 

In the same way, the difference mode is said to be in a 
two-mode squeezed state if either Ad+ < 1 or Ad 2 , > 1. 
On account of Eqs.fTUl) and (fTTj) along with the steady 
state solutions at the Appendix, the quadrature variances 
of take the form 

M 2 _ A 2 {1 - r? 2 ) + {2k + A + At])(2k + At? t Ay/1 - r? 2 ) 
±_ 2(2k + At))(k + Arj) 

(13) 

where rj = p^J — Pal which relates the initial probabil- 
ities for an atom to be in the upper and lower levels. 

According to Eq. ([9]) and (fT3|) . all quadrature vari- 
ances are independent of the parameter s which repre- 
sents the driving coherent fields. This shows that the 
driving coherent fields do not have any effect on the de- 
gree of squeezing of the two- mode light. This is due to 
the fact that the external coherent fields do not introduce 
additional coherence to the system which is believed to 
be the source of squeezing in three-level cascade lasers 
[IIMEEE!- Comparing Eqs. © and we easily 
to see that Ac\ = Ad 2 _ and Ac 2 _ = Ag?^ . It is found 
that exactly the same amount of squeezing exhibited in 
the minus quadrature of the sum mode and in the plus 
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FIG. 3: Plots of the quadrature variance Ac?, and Ad+ ver- 
sus 77 for k = 0.15 and for different values of the linear gain 
coefficient. 



quadrature for the difference mode. It is however worth 
mentioning that this statement is valid only at steady 
state. In the following we explore how the two mode 
squeezing can be optimized by varying the three param- 
eters in Eqs. © and (fl3"]l. In Fig. 2, we plot the Ac 2 _ 
and Ad+ versus cavity damping constant, n. As can be 
seen from this figure the variances increase almost lin- 
early with k, i.e, the squeezing decreases linearly with 
k. Obviously, k represents the property of one of the 
mirrors, port mirror, that builds the cavity. The larger 
the value of k the more the cavity modes susceptible to 
interact with the environment which essentially degrades 
the degree of squeezing. We also plotted, in Fig. 3, the 
variances of the squeezed quadratures versus the param- 
eter rj for different values of the linear gain coefficient. 
It is possible to see from this plots that the two-mode 
squeezing increases with the linear gain coefficient as pre- 
viously established [TJ, LUa, ll9(- Moreover, the value of 
, rj at which the maximum squeezing occurs decreases to 
zero as the linear gain coefficient increases. Therefore, 
two-mode squeezing can be optimized by choosing small 
valuesry (i.e, by initially prepaying slightly more atoms in 
the lower level) and large values of A/k. 



IV. ENTANGLEMENT BETWEEN THE TWO 
CAVITY MODES 

In this section we demonstrate a continuous-variable 
bipartite entanglement between the two cavity modes a 
and b. Entanglement is solely a property of quantum 
mechanics which is related to the inseparability of the 
combined density matrix of a system into the density 
matrices of the individual system. Inseparability criteria 
for continuous variable bipartite states have been pro- 
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posed by various authors 0, HH . In this work, we ap- 
ply entanglement criterion introduced by Duan et al. [9j 
which is sufficient and necessary condition for Guassian 
states and sufficient, in general to detect bipartite contin- 
uous variable entanglement. According to this criterion, 
a quantum state of a system is said to be entangled if the 
sum of the variances of the EPR-like quadrature opera- 
tors satisfy the inequality 



where 



and 



Au 2 + Aw 2 < 2{z 2 + 4), 



u = \z\x a + -X b 

z 



V = \z\p a p b 

Z 



(14) 



(15) 



(16) 



with x a — (a) + d)/v2, &b 



(St + b)/V2, p a = i(flt - 

a)/v2> and pb = i(w — b)j\f2 in which z is a non-zero 
real number. In the following, we choose z = — 1 so that 
the upper bound of the inequality Eq. (IT41 to be 2. This 
is not an optimal choice in general but it is sufficient to 
compare with other inequality. With z = — 1, the quadra- 
ture operators between which we want to demonstrate bi- 
partite entanglement are simply the squeezed quadrature 
operators for the sum and difference modes, u — c_ and 
v = d + . We already calculated the variances of these 
quadrature operators in the previous section. Thus on 
account of Eqs. © and (|T3|) , the sum of the variances of 
it and v becomes 

Au 2 + Av 2 = 2Ac 2 _ = 2Adl 



A 2 (l - rf) + (2k + A + At 1 ){2k + Aq - Ay/l-rf) 



(2k + Ar))(K + Arj) 



(17) 



We immediately notice that, this particular entangle- 
ment measure is directly related the two-mode squeezing, 
as previously reported elsewhere [HI, [l3j]. This direct 
relationship shows that whenever there is a two-mode 
squeezing in the system there will be entanglement in 
the system as well. It also follows that the degree of 
entanglement does not depend on the external driving 
coherent fields. This is attributed to the fact that the 
coherent fields do not introduce additional atomic coher- 
ence to the system, as the same is true for the case of 
squeezing. In Fig. 4, we show results for the sum of 
the squeezed quadrature variances indicating a clear vi- 
olation of the inequality and hence the entanglement be- 
tween the modes. It can easily be seen that the degree of 
entanglement increases with the rate at which the atoms 
are injected into the cavity, A. The entanglement how- 
ever disappears at two extreme value of rj, one at r\ = 
which corresponds to maximum injected atomic coher- 
1/2 and the other at 77 = 1 corresponds to 
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FIG. 4: Plots of Au 2 + Av 2 of the two-mode light in the cavity 
as steady state versus r] and for different values of the linear 
gain coefficient. 



V. MEAN PHOTON NUMBER OF THE 
CAVITY MODES 

In order to know how intense the produced light is, 
we calculate the total mean photon number of the cavity 
modes. In terms of the sum and difference operators 
of the cavity modes, the mean photon number can be 
written as 



(N) - (c f c) + (<?d). 



(18) 



Using the definitions of the sum and difference operators ( 
Eqs. (|6|) and (fT0]) h the mean photon number becomes 



(N) = (a f a) + (b r b). 



(19) 



Now, with the aid of Eqs. (|Al9j) and (|A20|) the mean 
photon number takes the form 



(N) = 



A(l -t))(2k + A + Arj) 



A 3 r)(l - i! 2 ) 



4(k 2 
Ae 2 [A 2 (\ - 



kAtj) 



4(k 2 + kAt])(2k + Arj) 



sjl^f) + 2(k 2 + kAt})} 



(k 2 + K,Ar\f 



(20) 



no atomic coherence, p£} = 



In Eq. ([20)1 the term that contain e represents the 
contribution from the external driving coherent fields to 
the total mean photon number. Graphically, the effect 
of the coherent fields on the mean photon number of the 
cavity modes is shown in Figs. 5 and 6. Fig. 5, indeed, 
clearly indicates that the mean photon number of the 
cavity modes increases with the amplitude of the coher- 
ent driving fields. In order to clearly see by what extent 
the coherent fields enhance the mean photon number over 
the laser system, we plotted in Fig. 6 the mean photon 
number versus the parameter rj in the absence and pres- 
ence of the coherent fields. It is quite interesting to note 
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FIG. 5: Plot of the mean photon number versus e and r\ for 
A/k = 100. 
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FIG. 6: Plots of the mean photon number versus r\ for and 
A/k = 100, in the absence of the coherent fields (solid) and 
in the presence of the coherent fields with e/k = 50 (dotted). 



from this figure that the coherent fields enhance the mean 
photon number over the laser system by several tens of 
thousands of mean photon numbers. More importantly, 
the increase in the mean photon number is observed in a 
region where the degree of two-mode squeezing and en- 
tanglement is significant making the system under con- 
sideration a viable source of intense squeezed as well as 
entangled light. 



driven by two external coherent fields at resonance with 
the cavity modes, where the atomic coherence is induced 
by initial superposition of top and bottom levels. We 
have demonstrated two-mode squeezing in the sum and 
difference fields with equal amount of squeezing in each 
field. In addition, we have shown that the two cavity 
modes are strongly entangled and the degree of entan- 
glement is directly related to the two-mode squeezing. 
We also found that the effect of the external coherent 
fields is to increase the mean photon numbers consider- 
ably in a region where the squeezing and entanglement 
are strong making the system under scrutiny a viable 
source of a bright macroscopic squeezed as well as en- 
tangled light. The degree of squeezing and entanglement 
can be optimized by choosing small values of the cavity 
damping constant, for large values of the linear gain co- 
efficient, and by initially preparing slightly more atoms 
in the lower level than in the upper level. 



Appendix A: Equation of evolution of cavity mode 
operators 

In this Appendix, we derive, applying the master equa- 
tion obtained in Sec. II, the equation of evolution for 
the cavity mode operators. We also calculate the steady 
state solutions of the resulting equations. The equation 
of evolution of an operator O in the Schrodinger picture 
is expressible as 



(Al) 



Using the master equation, Eq. (|5b| and Eq. ©, we 
obtained the following equations: 



-<a> = - T <a>-— <&t) +e , 



dV ' 2 W 2 



,&■>+£, 



(A2) 



(A3) 



-(a 2 ) = -Ma 2 ) ~ ApM{cft) + 2e(a), (A4) 



-(S 2 > = - M6 (6 2 } + ApW{tfb) + 2 e (6), (A5) 



-<fito) = -MaM - -^((&b) + (aW)) 
+ £ ((at) + (a))+Ap(° ), 



(A6) 



VI. CONCLUSION 

We have studied the steady state two-mode squeezing 
and entanglement in a two-mode three-level cascade laser 



|(6t6) = - M6 <6t6) + ^«a6) + <at6t)) 
+ e((b') + (b)), 



(A7) 



|(aSt) = -^t) + ^-((a 2 )-(^)) 
+ <(&) + (fo t », 



-(ab)=-»{ab) + -^({a)a)-{tfb)) 
+ \A P ^+e{{a) + (b)), 



in which 



A , = i[2 K + A(p(° c )-p( Q ))]. 



(A8) 



(A9) 



(A10) 



(aS+) = 



£ 2 (2/t-yl(l _^))(2/t + ^(l+>?)) 
(k 2 + kAt;) 2 



+ 

(aS) = 
+ 
+ 



e 2 A 2 y/T^f(2 - y/l^n 1 ) 



(k 2 + nArj) 2 



(A17) 



KA v / l-7 1 2 (2k + A(l + rj)) 

4(k 2 + kAti)(2k + Arj) 
e 2 (2n - A(l - t}))(2k + A(l + r,)) 
(k 2 + kAt)) 2 

e 2 A 2 ^T^f(2 - yT^ 2 ) 
(n 2 + kAt]) 2 



(A18) 



It proves to be convenient to introduce a parameter 
that relates the probability for an atom to be in the upper 
and lower levels, such that 



„(o) = l _A 

raa o ' 



JO) = l + V 

rcc n 



and 



(All) 



(A12) 



Next, we calculate the steady solutions by setting the 
time derivatives of the cavity mode variable to zero. By 
simultaneously solving the steady state solutions, we find 
the following: 



r . _ e[2« + A(l + 77)-yyr^2] 
W ~ n 2 + kAt) 



(A13) 



. t .. _ A(l -t])(2k + A + At]) 



(a T a) = 



4(k 2 + kAt]) 



(k + Ar!)A 2 (l-r, 2 ) 
4(2k + At))(k 2 + kAt)) 



4s 2 (2k + A(l + ??) - Ay/1 - rj 2 ) 
(2k - A(l - t]))(k 2 + kAt)) 

e 2 As/T^q 2 {2k + Ajl + rj)) 
(k 2 + kAi]) 2 



e 2 A 3 (l- V 2 )(2- y/T^y 2 ) 
(2k-A(1 -f]))(K 2 + kAt/) 2 ' 



(A19) 



(6> = 



e[2K-A(l-?7)-A v /l-r? 2 ] 
k 2 + kA?7 



(A14) 



_ 4e 2 (2 K + A(l+ V )-Ay/T^f) K A 2 (l-n 2 )) 
{a } (K 2 + K A v )(2n-A(l-r 1 )) {tfb) - - ' 1 1 



e 2 Ay/\ -rf(2n + A(l + ??)) 

(k 2 + kAt]) 2 
e 2 A 3 (l- v 2 )(2- y/T^f) 



4(2k + At])(k 2 + kAti) 



4e 2 (2n-A(l-r 1 ) + Ay/l-T 1 2 ) 



(2k + A(l + r]))(n 2 + kAt]) 

(A15) s 2 Ay/T^f(2k - A(l - r,)) 
(2. A(l V ))(n 2 + kA V ) 2 + {k2 - kAii)2 

+ f A \-™7f=?l (A20) 



p = 4e 2 (2 K - A(l - r,) + Vj - r? 2 ) 1 (2k + A(l + t?))^ 2 + K A V ) 2 

{ ' (k 2 + kAt))(2k + A(l + t?)) 

e 2 Ay/T^(2 K - A(l-rj)) 
+ (k 2 + nArj) 2 



-). £ Q !? V_j__^)_ (A16) It is worth mentioning that these solutions are obtained 

(2k + A(1+t]))(k 2 + kAj]) 2 ' provided that r, > 0. 
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